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Abstract 

A model of self-avoiding walk with suitable constraints on self-attraction is developed to describe 
the conformational behaviour of a single stranded short DNA molecule that form hairpin structure. 
Using exact enumeration method we calculate the properties associated with coil-hairpin transition. 
Our results are in qualitative agreement with the experiment. 
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The structural and dynamic behaviour of a single stranded DNA molecule that form a 
stem-and-loop (hairpin) structures in solutions is a subject of current interest (I], 0, f|. 
This is because DNA hairpin structure is known to participate in many biological functions, 
such as the regulation of gene expression |5[], DNA recombination |J and facilitation of 
mutogenic events J7J. The DNA hairpin-containing domains may provide potential binding 
sites for exogenous drugs and endogenous proteins ||. It can also be used as DNA biosenser 
{e.g. molecular beacons) ||. 

The single stranded DNA that form a stem-and-loop structures under a given set of solu- 
tion conditions are formed from synthetic oligonucleotides. For example, the oligonucleotide 
5'— CCCAA - (X) m - TTGGG— 3', where X is any one of the nucleotides and m its num- 
ber, is such a molecule which can easily be synthesized. Due to thermal fluctuations such a 
molecule may acquire different conformations. In a simplified description all of the configu- 
rations can be divided into two main states; the open and the closed one as shown in Fig. 
1. The closed hairpin structure is stabilized due to pairing of complimentary bases at the 
two arms of the molecule when hydrogen bonds are formed between them. The open state 
has high entropy due to the large number of configurations achieved by a single stranded 
DNA chain. The closed-to-open transition requires sufficiently large energy to unzip all the 
base pairs, whereas the closing transition requires the two arms of hairpin to come close to 
each other in space so that hydrogen bonding between the complimentary nucleotides can 
take place. Experimentally one finds the rate of closing depends strongly on the properties 
of the hairpin loop, such as the length and rigidity, whereas the rate of opening is relatively 
unaffected by these properties 0. 

In this note we describe a model and show that it explains some of the conformational 
behaviour associated with the formation of hairpin structure. 

We represent a linear polymer chain of N monomers by a self-avoiding walk of (N — 1) 
steps on a two dimensional square lattice. Visited lattice sites represent the monomers of 
the chain. The first n sites of the walk represent nucleotides, say, A or C and the last n sites 
the complimentary nucleotides, T or G. The remaining (N — 2n) = m sites represent the 
nucleotides of one type (either A or T or C or G) . The repulsion at short distance between 
monomers {i.e. excluded volume) interaction is taken into account by the condition of self 
avoidance. When the complimentary monomers in a walk occupy nearest sites in a particular 
direction (say, clock wise), they form a pair (known as base pair). The energy associated 
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FIG. 1: Schematic of the conformational fluctuations of DNA hairpin- loop. This molecular beacon 
can fluctuate between open and closed states in solution at ambient temperature. 

with this pair is e (e < 0). The base paring cannot take place if the nucleotides approach 
to each other in any other way. This condition has been imposed to take into account 
the fact that in a single strand DNA chain the pairing between complementary nucleotides 
takes place only when they approach each other directly without the sugar phosphate strand 
coming in between as shown in Fig. 2a. Fig. 2b shows the situation in which the sugar 
phosphate strand is in between and therefore pairing can not take place. 

All possible conformations of walks of N sites mapped by a self- avoiding walk are gen- 
erated using exact enumeration technique [|H], |ll|]. Since the time involved in enumerating 
these conformations increases as /i N where /i is the connectivity constant of the lattice, we 
have restricted ourselves to 12 < N < 22. 

The partition function of our interest is 

n 

Z N = Y,C N (p)(e^ k ° T y (1) 

p=0 

where p is the number of base pairs and CV(p) the total number of configurations corre- 
sponding to walk of (N — 1) steps with p number of pairs (base pairs). The Helmholtz free 
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FIG. 2: Schematic of configurations showing the (a) formation and (b) non-formation of base pairs, 
energy of the system is given as 



where ks is the Boltzmann constant and T the temperature of the system. From Eq. (2) 
one can calculate the specific heat and entropy using the following relations: 



In our calculation we have taken n = 5 and e = —0.08 eV. Dependence of the free energy 
(F) on temperature for iV = 22 is shown in Fig. 3a. 

From this figure it is clear that there is a "discontinuous" change in the free energy at 
T = 338°K. In Fig. 3b we plot the specific heat as a function of temperature for several 
values of N. This curve has a peak suggesting a transition from the hairpin structure to 
the coil (or open) state. For N = 22, the peak is at 337 ± 2°K. Since we are describing 
a finite system, the transition is not sharp; it takes place over an interval of temperature. 
This interval is obvious from the width of the peak in C v and also from Fig. 3c in which 



F = — 



k B T\nZ N (T) 



(2) 
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FIG. 3: (a) Plot of free energy (F) as a function of temperature (T) for TV = 22; (b) Specific heat 



we plot entropy (S) as a function of T for N = 22. The entropy is low at low temperature 
as the molecule is in a hairpin structural form whereas at high temperature the entropy is 
higher as the molecule in coil form assesses a large number of conformations. 

The partition function defined by Eq. (1) has six terms corresponding to six values of p 
from to 5 (p = corresponds to open state of the chain whereas 1 < p < 4 corresponds 
to partially closed state and p = 5 the final state of hairpin structure. We can find the 
probability of each of these states in the following way. 

P P {T) = f (4) 



where Z = J2 P =o Z p 

We compare P P (T) for different values of p in Fig. 4a. It shows that the dominating 
structure are either all closed or all open. 

The other quantity of interest is the fraction of open to closed structures at a given 
temperature. The fraction of all open structure can be derived from the relation 

X- = 7 ^° y (5) 
where Z is for p = and Z c = Ylp=i Z p an d the fraction of closed structures 



From Eq. (5) and (6) we get 

X- 



exp 

X+ 



1 IF -F C ) 



(7) 



where F Q = — k B T In Z and F c = — k B T hi Z c 

The fraction of all closed structure at a given temperature is found from the relation 

From this figure [Fig. 4b] one may note that while \- an d X+ intersects at T — 337±2°K 
in agreement with the coil-hairpin transition temperature, X- an d X+ intersects at T — 
366 ± 2°K. As shown in Fig. 4a the probability of partially connected conformations are 
small and confined to narrow temperature region. One may therefore conclude that the 
conformational properties of this kind of single DNA strand is primarily controlled by either 
all open or all closed states. 
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FIG. 4: (a) Probability of opening (p = 0), closed (p = 5) and partially closed (p = 1, 2, 3,4) states 
as a function of temperature (T); (b) (i) Fraction of all open structure (x~), (ii) fraction of closed 
structure (x+) an d (iii) fraction of all closed structure (x'+) as a function of temperature (T). 

By varying the value of e, one can study the effect of base pair interaction on coil-hairpin 
transition. The change in e value from -0.08 eV to -0.05 eV leads to the decrease in transition 
temperature by 100°K. 

In conclusion, we have studied a simple model for DNA hairpin loop structure and their 
relative stabilities. We use lattice models of polymers along with exact enumeration tech- 
nique to study the stability of different conformations of DNA hairpin loop and the corre- 
sponding hairpin-coil transition. Our results are in qualitatively agreement with the reported 
experimental work by Bonnet et al [|1|]. 
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